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We have started a program to compute the electromagnetic form factors of mesons. We discuss the techniques 
used to compute the pion form factor and present preliminary results computed with domain wall valence fermions 
on MILC asqtad lattices, as well as Wilson fermions on quenched lattices. These methods can easily be extended 
to p — *• 77r transition form factors. 



1. INTRODUCTION 

The pion electromagnetic form factor is of- 
ten considered a good observable for studying 
the relative accuracy of perturbative and non- 
perturbative descriptions of QCD as the energy 
scale varies. It is hoped that because the pion is 
the lightest and simplest hadron that perturba- 
tive descriptions will remain accurate at lower en- 
ergy scales than predictions for heavier and more 
complicated hadrons, e.g. the nucleon. 

At very low energies, the charged pion moves 
moves in an electromagnetic field like a point par- 
ticle with one unit of electric charge. This allows 
us to fix the overall normalization of the form fac- 
tor 

F X (Q 2 ) = 1 as Q 2 ^0. (1) 

The experimentally observed behavior of the form 
factor at small momentum transfer is accurately 
described by the vector meson dominance (VMD) 
hypothesis |T12I3| 

Ftt(Q 2 ) ~ - * , 2 for Q 2 «m 2 (2) 

The current experimental situation is presented 
in Fig. ^ taken from 0]. 
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sources at ORNL, under the auspices of the U.S. DoE's 
SciDAC initiative. 
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Figure 1. Summary of experimental data and 
phcnomenological predictions for the pion elec- 
tromagnetic form factor, taken from 0]. Note 
that some of the points show projected errors for 
proposed experiments. 

What is surprising is that VMD with only the 
lightest m p resonance appears to describe all the 
existing data, even up to scales of Q 2 > lGeV 2 
where one might have hoped the perturbative 
QCD (pQCD) would provide accurate predic- 
tions. At very high momentum transfer, we ex- 
pect the perturbative description will be correct 
as computed in |5I6I7| 

F,m ^ a fP fl as Q^oo (3) 

While the data may indicate an onset of the cor- 
rect scale dependence around 2 GeV 2 , the numer- 



ical value is 100% larger than the pQCD asymp- 
totic prediction. 

This situation presents many questions. At 
what scale does the form factor scale with Q 2 as 
predicted by pQCD? Even if we assume proper 
scaling, the numerical values do not agree with 
pQCD. How rapidly will the data approach the 
pQCD prediction and at what scale will pQCD 
finally agree with the data? These are questions 
which Lattice QCD calculations are ideally suited 
to address, provided we can get reliable results for 
momentum transfer on the order of a few to sev- 
eral GeV 2 . 

Early lattice calculations validated the vector 
meson dominance hypothesis at low Q 2 |8I9I . Re- 
cent lattice results |10lll| , including some of our 
own preliminary results |12| . have somewhat ex- 
tended the range of momentum transfer, up to 
2 GeV 2 , and the results remain consistent with 
VMD and the experimental data. 

2. LATTICE COMPUTATION OF F^{Q 2 ) 

The electromagnetic form factor is obtained 
in lattice QCD simulations by placing a charged 
pion creation operator at Euclidean time ti, a 
charged pion annihilation operator at tf and a 
vector current insertion at t as shown in Fig. [2] 
A standard quark propagator calculation pro- 
vides the two propagator lines that originate 
from ti. The remaining quark propagator, orig- 
inating from tf is obtained via the sequential 
source method: (1) completely specify the quan- 
tum numbers, including momentum pf, of the 
annihilation operator to be placed at tf and (2) 
contract the propagator from ti to tf to the an- 
nihilation operator and use that product as the 
source vector of a second, sequential propagator 
inversion. The resulting sequential propagator 
appears as the thick line in Fig. extending from 
ti to t via tf. Given these two propagators, the 
diagram can be computed for all possible values 
of insertion position t and source and insertion 
momenta, pi and q, consistent with overall mo- 
mentum conservation: q = Pf — pi- Furthermore, 
with the same set of propagators, any current can 
be inserted at t and any meson creation oper- 
ator can be contracted at ti. So, the diagram 



relevant to determining the form factor for the 
transition p + — > jir + can be computed without 
further quark propagator calculations. By apply- 
ing the sequential source method at the sink, the 
trade-off is that the entire set of sequential prop- 
agators must be recomputed each time new quan- 
tum numbers are needed at the sink, particularly 
Pf. 




t.,p. 



Figure 2. The quark propagators used to com- 
pute the pion form factor. 

We can extract the pion energies E 7r (p) using 
standard lattice techniques of fitting pion corre- 
lation functions from which we can compute the 
momentum transfer 



- Q 2 = [E„(p f ) - E„(pi)} 2 - (p f - pif 



(4) 



which should be non-positive if the pion spectral 
function is well-behaved. Since the largest Q 2 oc- 
cur in the Breit frame, Pf — ~Pi, it is important 
to choose a non-zero pf to achieve large momen- 
tum transfer. 

The form factor, F(Q 2 ), is defined by 

|^(0)| 7T(p i )> continuum = (5) 

Z V (n(p f ) |^(0)| 7r(p-)> = F(Q 2 )( Pl + P/ ) M 

where Vn(x) is the chosen vector current. The 
three-point correlation function depicted in Fig. [21 
is given by 



^n^( t i,t,tf ) p i ,Pf) 



(6) 



a 9 e ~i(x f -x)-pf x e -i(x-Zi)-pi 



x (o|^(i/)^(s)^(n) 



where 4>\{x) and 4>b(x) are creation and annihila- 
tion operators with pion quantum numbers. The 



A and B indicate that different operators may be 
used at the source and sink, i.e. smeared source 
and point sink or pseudoscalar source and axial 
vector sink. 

Inserting complete sets of hadron states and 
requiring tj -c t <C tf, gives 



(7) 



a 3 (0|<M*)k(P»> 
x M£f) \V^{x)\ir(pi)) 



2E 7! {p l ) 

Similarly for the two-point correlator and ti -C tf 
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(U,t fl p)^a 3 (0\<l> B (xf)\nffi) 

e -(tf-ti)E«(p) 



(8) 



For pion operators, we use either the pseu- 
doscalar density, c/)^(x) — ip{x)^f^{x), or the 
temporal component of the axial vector current, 
4>( 2 \x) — ^(x)7574^i(x), at a single point, de- 
noted "L" for local, or smeared, denoted "S," us- 
ing gauge invariant Gaussian smearing 
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(9) 



Note that we have suppressed the flavor struc- 
ture. For these operators, we get the following 
amplitudes 
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(13) 



We anticipate that Z L (p) oc E n (p) based on the 
Lorentz structure of the operator. All of the Z's 
should have corrections of 0(p 2 ). 

From this discussion, we can see that there are 
two ways to determine the form factor F W (Q 2 ). 
The first method, which we will call the fitting 
method involves a fit of the relevant two and 
three-point functions to simultaneously extract 



Table 1 

Simulation details for quenched Wilson fermion 
calculations at a -1 » 2 GeV 

k volume aniir am p 

0.1480 16 3 x 32 0lS73 0.712 

0.1520 16 3 x 32 0.477 0.549 

0.1540 16 3 x 32 0.364 0.468 

' 0.1555 24 3 x 32 0.259 0.398 

0.1563 24 3 x 32 0.179 0.358 

0.1566 24 3 x 32 0.145 0.343 

0.1566 32 3 x 48 0.145 0.343 

the form factor, the energies E n {p) and the am- 
plitudes Z(p) in a single covariant, jackknifed fit. 

Another method, which we will call the ratio 
method starts by determining the energies E 7r (p) 
and then constructing the following ratio which 
is independent of Z^\ Zs{\p\) and all Euclidean 
time exponentials: 



F(Q 2 ,t) = 



r£? ff (u,t,t f ,p h p f )T^ (u,t, pf ) 



2Z v E v (p f ) 



(14) 



E^pi + E w (ff) t 

where the indices A, B and C can be either L 
(local) or S (smeared). For a conserved current, 
Zy = 1. As part of our program, we expect to 
determine the relative merits of each extraction 
method. 

3. SIMULATION DETAILS 

Our first calculations were done on quenched 
configurations generated with the Wilson gauge 
action at (3 = 6.0 (a -1 « 2 GeV). The propa- 
gators were computed using the unimproved Wil- 
son fermion action with Dirichlet boundary con- 
ditions. For Wilson fermions at these lattice 
spacings, the exceptional configuration problem is 
rather mild particularly when compared to non- 
perturbatively improved Clover action. This en- 
abled us to reach pion masses of 300 MeV with- 
out observing any exceptional configurations, see 
Tab.n wnereas Clover simulations are limited to 
pions of roughly 500 MeV |10|. Of course, pion 
masses in quenched domain wall fermion calcula- 
tions arc limited only by finite volume effects and 
available computing power. Yet the best results 



so far for pion form factors are 390 MeV pions at 
a -1 w 1.3 GeV lattice spacing Furthermore, 
at this workshop we learned that Wilson fermion 
results may be automatically 0(a) improved with 
just double the effort |13j . which is an option we 
have under consideration for the future. 



Table 2 

Simulation details for domain wall fermion calcu- 
lations on 272 dynamical MILC asqtad lattices 
at a- 1 w 1.5 GeV. 
am u d am s amawf Volume am n 



0.01 



0.05 0.01 



20 ]i x 64 0200" 



For our unqucnchcd results, we used valence 
domain wall fermions with a domain wall height 
of too = 1-7 and L s = 16 on MILC N f = 2 + 1 
and Nf = 3 lattices after HYP blocking [H] . 
Dirichlet boundary conditions were imposed 32 
timeslices apart for the domain wall propagator 
calculations. Detailed results of several quantities 
computed on these MILC configurations were pre- 
sented at this workshop |15|. Preliminary results 
on related hadronic observables using many of 
the same domain wall propagators have recently 
been presented ^B] and were also discussed at this 
workshop. 

4. PRELIMINARY RESULTS 

Representative plots of the pion form factor 
and it's dependence on the timeslice of the current 
insertion ^(Q 2 , t) are shown in Fig. EI The form 
factors in these plots were computed using the ra- 
tio method. For the heavier pion masses (smaller 
values of k) there are nice plateaus which indicate 
reasonable signal to noise for all momenta consid- 
ered and little contamination from excited states. 
As the pion masses get lighter, the signal to noise 
decreases and the plateaus become less convincing 
which suggest that the source and sink may not 
be sufficiently separated and hence excited states 
may be contaminating the ratio. Another possi- 
bility is that the ratio may not provide an optimal 
determination of the form factor at lighter pion 
masses. A detailed study is currently in progress. 

Preliminary results of the quenched Wilson 
form factor computed by the fitting method are 



plotted in Fig. 21 Also shown in the figure are ex- 
perimental data points and a curve showing the 
prediction of vector meson dominance using the 
observed value for the p meson mass. While the 
data tend in the correct direction with decreas- 
ing pion mass, the reader may notice that the 
form factor for 300 MeV pions already lies below 
the physical curve. This suggests that fitting the 
data to extract a lattice determination of the vec- 
tor meson mass would underestimate the physi- 
cal value. In fact, this is exactly what one should 
expect since it is known that 0(a) scaling viola- 
tions to the vector meson mass computed with 
Wilson fermions tend to underestimate the mass 
by roughly 20% J7|, the same amount needed to 
move the form factor points above the continuum 
curve. Hence, the proposal of Frezzotti and Rossi 
[T3*j becomes all the more appealing as a straight- 
forward way to remove what is presumably an 
0(a) lattice artifact. 

Because we would like to compute the pion 
form factor at large momentum transfer, we have 
spent a substantial amount of effort on our do- 
main wall data set in extracting the pion energies 
and amplitudes at relatively large momenta. In 
the continuum limit, the pion dispersion relation 
should follow the continuum one 



(15) 



Another possibility is that the pion dispersion re- 
lation will follow the dispersion relation of a free 
lattice boson 



— sinh 2 



(16) 
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Both relations agree in the small momentum 
limit. 

In Fig. 03 we have plotted against each dis- 
persion relation and we see that both dispersion 
relations provide a reasonable representation of 
the data, although there may be a slight flat- 
tening of the data against the continuum curve 
at higher momenta. These results suggest that 
directly fitting all the data to either dispersion 
relation, thereby reducing the number of fit pa- 
rameters needed to extract the form factor, may 
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Figure 3. Pion form factor data versus timestep, for (a) k=0.1520, (b) k=0.1540 and (c) re=0.1563. 
Numerical values of Q 2 in GeV are shown in the legends. 
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Figure 4. Results for the pion form factor as a function of Q 2 for each of the available k values. Experi- 
mental measurements [T%] and the vector meson dominance hypothesis are also shown. 



improve the relative signal to noise of the remain- 
ing parameters. This may help dramatically in 
the ratio method, where the only fit parameters 
are the form factor and the energies. 

In the fitting method, one must reliably esti- 
mate not only the energies, but the amplitudes, 
at high momenta. In Fig. [5] we present the four 
amplitudes we estimate from the four two-point 
correlators we measure: smeared-smeared and 
smeared-local for both pseudoscalar-pseudoscalar 
and axial-axial operators. In the fitting proce- 
dure, all four correlators are constrained to have 
the same energy. From the figure, we can see that 



our expectations of Z L oc const and Z L (p) oc 
Ek{p) are consistent with the data. We can also 
see from Z^ that the smeared pseudoscalar op- 
erator has a strong overlap with the zero momen- 
tum pion but that the overlap diminishes rapidly 
with increasing momenta. From Z s we see that 
the smeared axial-vector operator has good over- 
lap over a wide range of non-zero momenta, with 
maximal overlap around ap ~ 0.5. Thus, given 
the variety of sources and their relative strengths 
of overlapping with pions at different momenta 
was key to fitting pion energies up to such high 
momenta. 
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Figure 5. Pion dispersion relation vs. continuum (left) and lattice (right) expectations. 
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Figure 6. Results for the four pion source amplitudes Z^(p) used in this study. 



5. CONCLUSIONS 

From our preliminary quenched Wilson form 
factor results, we find that both the ratio method 
and the fitting method are useful tools for com- 
puting the pion form factor. Each method has 
different systematic errors, so the extent to which 
both agree should give confidence that the sys- 
tematic errors are small and well understood. 

From our preliminary dynamical domain wall 
spectrum results, we recognize the importance of 
using a large basis of pion operators so that at 
least one will have reasonable overlap with the 
momenta under consideration. The local axial 
vector operator is particularly useful in fitting 
higher momentum states since it's overlap with 
a given state increases as the energy increases. 

In the near future, when the pion form factor 
analysis is complete, it will be a straightforward 
extension to our existing analysis framework to 
compute the p + — > -fn + transition form factors 
using the propagators and sequential propagators 
we have already computed. 
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